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Recent results on the hadroproduction of charm and beauty are discussed in the 
context of the current theoretical and experimental status. The areas covered 
include production of open charm and beauty particles and charm and beauty 
quarkonia. Experimental results from both fixed target and colliding beam exper- 
iments are presented and compared to theoretical predictions. Predictions based 
on QCD perturbation theory are generally successful in explaining the shapes of 
differential cross sections but underestimate the scale of the cross section by factors 
of 2-3. An exception is production of the vector meson quarkonia which cannot be 
satisfactorily explained by perturbative processes alone. 



1 Introduction 

This review emphasizes recent progress in the understanding of the hadropro- 
duction of charm and beauty with an attempt to explain the new developments 
in the context of the fundamental questions of interest. Experimental and the- 
oretical work is presented which best illustrates the subject. In cases where 
two or more illustrations of the same topic are available, usually one is cho- 
sen for presentation here and the others cited as references. In making these 
choices, some preference is given to results published in the regular journals 
over preliminary results. However, important preliminary results are also pre- 
sented here. As there are many recent publications in this field, the author 
sincerely hopes that the fraction of the available results presented here serves 
to stimulate interest in the whole. 



2 Theory 



At high energy, heavy quark production is dominated by gluon collisions. For 
example, in a pp collision, illustrated in Figure 0, the valence quarks are spec- 
tators while the hard interaction producing the heavy quark pair is via gluon 
fusion. Many important characteristics of the heavy quark production proper- 
ties are evident from consideration of the gluon-gluon collision in its center of 
mass frame as illustrated in Figure |^. At high energy, the g — g CM energy is 
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given by 

S = X1X2S (1) 

where Xi are the momentum fractions of the gluons and s is the hadron-hadron 
CM energy. Since s must be above threshold for a pair of heavy quarks, 
the cross section clearly depends on the gluon distribution functions of the 
colliding hadrons. Thus, pions, kaons and protons incident on the same target 
will produce different heavy quark cross sections. The g — g CM frame is in 
general moving with respect to the hadron CM frame, and the heavy quarks 
can have substantial longitudinal boosts in the hadron CM frame. The more 
unequal the gluon momentum fractions, the larger will be this longitudinal 
boost. Since the gluon distribution functions peak at small x, the heavy quark 
longitudinal momentum distributions (in the hadron-hadron CM) will peak at 
zero and fall rapidly with increasing magnitude of the longitudinal momentum. 
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Figure 1: Feynman diagram illustrating a pp collision producing a heavy quark pair 



The HQ production cross section is given bylJ, 

a = J dxidx 2 a{s,q 2 ,MQ)g l {x u ql)g 2 {x2,ql) (2) 

where gdxi,qf) are the gluon distribution functions for momentum fraction 
Xi and momentum transfer qf, and a is the short distance cross section. The 
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Figure 2: The gluon gluon collision in its center of mass frame. 



strong coupling constant appears in a and has the form, 
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Perturbation theory can be used to evaluate a when the coupling constant is 
less than one, as occurs when q 2 , which is of order Mq , is greater than &q CD . 
This is the case for the charm and beauty quarks. 

However, there are a number of theoretical problems with the formalism 
of equation 2. In existing calculations, the gluon distribution functions do not 
include heavy quarks in their evolution. Second, there are terms in the short 
distance cross section with coupling constants that become larger than one in 
certain kinematic regions. For example, factors of the form, 



become large for s > q 2 . For existing calculations, the theoretical problems 
are manifested by large sensitivity to changes in the renormalization and fac- 
torization scales, as discussed later. 

If the produced heavy quarks are close enough in momentum, they can 
form a bound state: the quarkonia states of charm or beauty. This formation 
process has a small probability and hence the fraction of heavy quark produc- 
tion into quarkonia is small, of order 1 part in 1000. The bulk of heavy quark 
production goes into open charm or beauty states. These cross sections de- 
pend on the formation of heavy quarks into particles, a process often referred 
to as hadronization or fragmentation. Several sub-processes can contribute to 
this hadronization, as illustrated in Figure [| The heavy quarks can interact 
with light quark/ anti-quark pairs arising from virtual gluons. Also, the heavy 
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quarks can interact with light valence quarks from the colliding hadrons which 
are spectators to the hard interaction which produced the heavy quark pair. 
The hadronization process can involve many small q 2 interactions and hence 
is intrinsically non-perturbative. Phenomenological models, such as the Lund 
string model, have been developed to describe the hadronization process. 

Because the amount of heavy quark fragmentation into a given charm 
species cannot be reliably treated theoretically, comparison of the predicted 
total heavy quark cross section with experiment must include measurement of 
as many species (mesons and baryons) as possible. Comparison of measured 
differential cross sections for particles with theoretical predictions for quarks 
can reveal the effect of fragmentation. 

As a preliminary remark, we note that rapidity, y, and Feynman-x, xf, 
are commonly used as variables to measure the longitudinal momentum of the 
produced heavy particle. Rapidity is the standard variable of special relativity. 
Feynman-x is defined, in the CM frame of the colliding hadrons, as the ratio 
of the 3-momentum of the heavy particle (or quark) to the maximum possible 
3- momentum. Hence, the range of Xp is from -1 to 1. 

3 Open Charm 

There are new measurements, of the production of open charm states from 
Fermilab Experiment E769B'tI This experiment utilized a 250 GeV charged 
secondary beam with the identified beam particles 7r + , 7r~, K + , K~ , and p. 
The charm particles D + , D°, D* + , Dg and A c are measured using fully recon- 
structed decays to all charged final states. Examples of mass plots are shown 
in Figure ^. The signal to background ratios and signal sizes are sufficient to 
measure the total and differential cross sections for most of the charm particles 
and beam types mentioned above. In some cases, only upper limits on cross 
sections are possible. 

E769 is a fixed target experiment utilizing a forward spectrometer. The 
acceptance for reconstructed charm particle decays is typically over the range 
— .1 < xf < -9. Charm particle cross sections are reported for the xp > 0, a 
range convenient for comparison with other experiments. Also convenient for 
comparing different experiments and theory is the sum of the inclusive cross 
sections, 

a(D+) + a(D-) + a(D°) + <r(E°) + <r(Df) + <r(Dj) (5) 

for xf > 0. This sum provides a partial measure of the total charm cross 
section for producing a pair of charm and anti-charm quarks. 
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Figure 3: The hadronization (also called fragmentation) process. Light quark pairs 
can form "strings" between the heavy quarks and between a heavy quark and a 

spectator quark. 
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Figure 4: (from ref. 2) Invariant mass distributions from Fermilab E769. The decay modes 
(including charge conjugates) and number of signal events with statistical errors are (a) 
D+ -» K-7T+7T+, 994 ± 47, (b) D° A"-7r+, 847 ± 55, (c) D+ -> E"+^-tt+ (combined </> 
7r and K*K modes), 100 ± 15, (d) — » pK~7r+, 35 ± 9. The signals shown are combined 

from all beam types. 
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Measurements of summed cross sections in equation for 7r~ -nucleon 
interactions from E769 and other experiments are plotted versus beam energy 
in Figure ^. Also shown are theoretical predictions based on perturbative 
QCD for different values of the renormalization scale, [1r. The predictions 
are for twice the cross section for producing a charm quark pair for which at 
least one of the pair has xp > 0. These predictions are generated using the 
computer program of Mangano, Nason and Ridolfi. It is clear from Figure || 
that the measured shape of the energy dependence is in reasonable agreement 
with theory and that the scale of the measured cross sections are higher than 
the central value of the prediction. Keep in mind that the summed cross 
section does not include charm baryons, which would increase the discrepancy 
if included. As can be seen from the curves for different values of /ir, the 
prediction has significant variation. This is also the case for variation of the 
factorization scale. 

For p nucleon interactions, the comparison of the summed cross section 
with theory is shown in Figure [| The comparison shows the same qualitative 
behavior as for the 7r~-nucleon case. 

There are no theoretical predictions for the kaon- nucleon cross sections 
since no parton distribution functions are available for this case. Cross sections 
for five different charm species and for pion, kaon and proton beams, measured 
by E769, are compared in Table [l]. Generally, for a given charm species, 
the cross sections are similar for the different beams. One notable exception 
appears to be the relatively large cross section for A c production with protons 
as compared to pions and kaons. This effect may indicate that an initial state 
light baryon enhances the production of a charm baryon. An analogous test 
is whether an initial state strange particle (kaon) enhances the production 
of charmed strange mesons (Df ). From the table, one sees no indication of 
this effect: the kaon induced D s cross section is consistent with the D s cross 
sections for pions and protons. 

The relative size of the A c and D s cross sections to those of the light 
charmed mesons is important for measuring the total charm cross section. 
The data of Table |l| from E769 indicate that the cross sections for A c and D s 
production are comparable to those of D + and D°. This is also supported by 
earlier data from CERN experiment NA32 (ACCMORp which are compared 
to those from E769 in Table § Thus, it appears that all charm species have 
comparable cross sections. The summed cross sections for all charm particles 
is at least a factor of three more than the central NLO QCD prediction. 

Another new result from E769 is the measurement of differential charm 
cross sections in xf and Pt for pion, kaon and proton beams S. As we will 
see, the cross sections are affected by the parton distributions in the different 
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Figure 5: (from ref. 2) Measurements of the summed cross section in equation (|5|) 
for 7r~ nucleon interactions from E769 and other experiments plotted versus beam 
energy. The NLO QCD predictions are for the sum of the charm quark and anti-charm 
quark cross section for xf > 0. The upper and lower curves result from variation of 
the renormalization scale from m c down to m c /2 and up to 2m c , respectively. 
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Figure 6: (from ref. 2) Measurements of the summed cross section in equation (g) for 
p nucleon interactions from E769 and other experiments plotted versus beam energy. 
The NLO QCD predictions are for the sum of the charm quark and anti- charm quark 
cross section for xf > 0. The upper and lower curves result from variation of the 
renormalization scale from m c down to m c /2 and up to 2m c , respectively. 
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Table 1: Charm particle (including anti particle) cross sections in /xb/nucleon for xp > 

measured by E769. 







A* 


P 


D+ 


3.2 ±0.2 ±0.2 


3.0 ±0.3 ±0.2 


3.2 ±0.4 ±0.3 


D° 


7.2 ±0.5 ±0.4 


7.2 ±1.0 ±0.4 


5.6 ±1.3 ±0.5 


D s 


2.0 ±0.4 ±0.2 


3.0 ±0.8 ±0.3 


> 0.5, < 2.5 


D*+ 


2.8 ±0.3 ±0.2 


1.7 ±0.5 ±0.1 


1.8 ±0.6 ±0.1 


A c 


3.3 ±1.1 ±0.5 


1.8 ±0.6 ±0.1 


> 5.0, < 21.2 



Table 2: Comparison of cross sections (including anti- particle) in ^tb/nucleon for xp > 
measured by NA32 (ACCMOR, 230 GeV tt~ beam) and E769 (250 GeV 7r± beam ). 





NA32 


E769 


D s 
A c 


1.4 ±0.2 ±0.2 
4.1 ±0.5 ±0.7 


2.0 ±0.4 ±0.2 
3.3 ± 1.1 ±0.5 



beam particles. 

For the E769 measurements, in each bin of the variable of interest (xf or 
Pt)i the cross section is summed as in equation (^J). (The summation increases 
the statistical power of the measurement.) The measurements for da/dxF are 
shown in Figure ^. Also shown are the NLO QCD predictions for charm quark 
production from pion and proton beams. (The predictions are normalized to 
best fit the data; the shape of the prediction has no constraints from the data.) 
It is clear that the expected shape for the proton beam is steeper in xf than 
for the pion beam. Indeed, the tt and p data are not consistent with having 
the same shape (C.L. > 99%). Furthermore, the tt and p data are well fit by 
the predicted shapes (C.L. > 50%). Thus, it appears that the Xf dependence 
shows the expected effect of the beam parton distribution. 

Although there is no prediction with a kaon beam, the data for pion and 
kaon beams are consistent with having the same shape (C.L. > 95%). This 
indicates that the pion and kaon have similar parton distributions. 

The differential cross sections in are shown in Figure |^. The pion 
beam data are consistent with the NLO QCD prediction for a pion beam, but 
inconsistent with the prediction for a proton beam. The difference in the Pt 
distributions for pions and protons is evident, although not as distinct as that 
for the xf distributions. 

One of the puzzles of charm particle production is that shapes of the dif- 
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Figure 7: (from rcf. 3) E769 measurements of da/d,XF for charmed mesons produced 
with 7r, K and p beams. Only statistical errors are shown. Arrows indicate 90% 
confidence level upper limits. The predictions are for charm quarks, based on NLO 

QCD. 




Pt 2 (GeV 2 ) 



Figure 8: (from ref. 3) E769 measurements of da/dP^ for charmed mesons with 
if > for charmed mesons produced with n, K and p beams. The predictions are 
for charm quarks, based on NLO QCD. 
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ferential cross sections in xf and Pt are so close to the NLO QCD predictions 
for charm quarks. This is certainly not the case in e + e~ production, where the 
effect of fragmentation is substantial. Indeed, if the Peterson fragmentation 
function is used to modify the NLO QCD result for charm quark hadroproduc- 
tion, the shape of the prediction falls much too rapidly compared to the data 
as shown in Figure |^. Perhaps loss of charm quark momentum from fragmen- 
tation is counterbalanced by an acceleration from interaction with spectator 
quarks. 




Figure 9: The solid curve is a fit to an early measurement of the D meson cross 
section versus xf from E769 with a pion beam. The dashed curve is the NLO QCD 
prediction for charm quark production. The dotted curve results when the Peterson 
fragmentation function is convoluted with the NLO QCD prediction for quarks. 



4 Open Beauty 

In contrast to measurements of the charm cross section, most measurements 
to date of the beauty cross section rely on partial reconstruction techniques, 
as did charm measurements of a decade ago. Many of these earlier charm 
measurements suffered from systematic errors due to the assumed models of 
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production and decay dynamics. This history provides a lesson for the current 
measurements of beauty production. For the measurements based on partial 
reconstruction techniques, systematic errors in the beauty quark pair cross 
section are difficult to avoid from model dependent assumptions about the 
production and decay dynamics. 

4-1 Fixed Target Beauty Experiments 

A variety of partial reconstruction techniques have been employed to measure 
the bb cross section in fixed target experiments. The experiments and tech- 
niques employed are listed in Table ||. The measurements, which are plotted 
in Figure [Io[ have large uncertainties, but generally lie at or above the central 
values of the QCD prediction. 



Table 3: Fixed target experiments measuring the beauty pair cross section. 



Experiment Name 


Technique 


CERN NAIO 
FNAL E653 
FNAL E672/E706 
FNAL E789 
CERN WA92 


tri muon 
nuclear emulsion 

inclusive J/ip with detached vertex 
inclusive J ftp with detached vertex 
silicon active target 



4-2 Measurements at the Tevatron 

At the Fermilab Tevatron collider, a well established technique to measure the 
beauty cross section is via th£ inclusive cross section for single high P t lep- 
tons. A recent measurement by the DO experiment, using muons, is shown 
in Figure [ll] as well as earlier measurements from CDF. Plotted is the cross 
section for a b quark as a function of the minimum pt of the b quark. The 
plotted cross section is derived under the assumption of the ISAJET model 
relating the pr of the lepton to that of the parent b quark. The model takes 
into account the momentum fragmentation of the b quark, the speciation into 
hadronic particles and the decay kinematics of those hadrons. Since few of 
these processes have been measured, systematic errors are inherent. The b 
quark cross sections extracted from both experiments (using the same produc- 
tion model) are consistent with theory, but fall along the upper range of the 
theoretical uncertainty. 
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Figure 10: Fixed target B cross section measurements. The solid lines are the 
predictions for a proton beam and the dotted lines are the predictions for a pion 
beam. The predictions are based on NLO QCD (ref. 6). Theoretical curves are 
shown for central parameter values and variations of Aqcd = 260 ± WOMeV and 

fi = fio, fJ-o/2, and 2/i . 
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Figure 11: (from ref. 7) Measurements of the cross section for a b quark derived from 
the spectrum of single leptons under the assumption of the ISAJET model. The solid 
and dotted lines show the prediction and uncertainty based on NLO QCD. For the 
prediction, the central values of parameters are Aqcd = 140 MeV and nit = 4.75 
GeV with fj,F = fJ-R = fJ.o = \J ml + {p b T ) 2 - The curves showing the theory error arc 
based on variations of AqcdI^q MeV, fM> —> fJ-o/2 and fj, — > 2/i . 
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Beauty cross section measurements based on exclusive decay modes of 
specific particles avoid some of the model dependence of the—measurements 
derived from inclusive lepton production. The CDF experimentO has measured 
the cross sections for the production of B + and B° through the decay modes 
B + — > tpK + and B° — » tpK*° with K*° — > if+7r~. Invariant mass plots for 
these decays, for different momentum ranges of the parent B, are shown in 
Figure |l2|. Analysis cuts include those based on the J/tp and K* masses and 
proper lifetime greater than 100/zm. 

The signals for the two decay modes in Figure [l^ are combined to measure 
an "average B meson" differential cross section, shown in Figure [13[ The 
measured differential cross section is consistent with the upper range of the 
theoretical prediction. Thus, this result shows the the same qualitative effect 
as the cross section extracted from the inclusive lepton method. Integrated 
over the range Pt > 6 GeV and for —1 < y < 1, the average B meson cross 
section is 2.4 ± 0.3 ± 0.4 /jb. 

Yet another type of measurement is available as a check on the inclusive 
lepton and exclusive decay methods discussed above. This (third) method is 
the measurement of the associated bb pair where the muon is tagged from one 
b quark and the jet is tagged from the other. This process has a distinct signa- 
ture in the plane transverse to the line of collision. The b and b quarks emerge 
nearly back to back, the deviation from collinearity being due to the transverse 
momentum of the initial state partons, for example, the gluon pair. When the 
b quark fragments, the muon tends to retain the direction of the b quark. Thus, 
the muon and jet from the accompanying b quark have a relative angular dis- 
tribution, Scj), in the transverse plane that peaks_at ir radians. The differential 
cross section for 5<p from the CDF experiment □ is shown in Figure The 
shape of the cross section is in reasonable agreement with the QCD prediction 
shown, but the data lie above the QCD prediction in absolute value. Thus it 
appears from a variety of measurements that the available QCD predictions 
underestimate the absolute bb cross section while being in good agreement with 
the shape of differential cross sections. 

Measurement of the total beauty cross section requires measurements for 
both beauty mesons and baryons. Information o p ft baryon production is just 
now becoming available. The CDF collaborationE3 has studied the decay chain 

A b -> J/tp A, J/tp -» A -> pn~ . (6) 

The charged tracks are required to have pt > 0.4 GeV. Vertex requirements 
are placed on the A, J/tp, and Ab. Cuts are also placed on the reconstructed 
masses of the A and J/tp. A signal of 8 ± 4 events is visible in the invariant 
mass plot for Ab shown in Figure EBI, on the left, for muon tracks measured 
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Figure 12: (from ref. 8) Invariant mass plots from the CDF experiment for (left 
column) B+ -> and (right column) B° -> tpK*° with if* -> if + 7r~. The 

Pt ranges are (a,e) 6-9 GeV/c, (b,f) 9-12 GeV/c, (c,g) 12-15 GcV/c and (d,h) > 15 
GeV/c. The signal in (a) is 51 ± 10 events and in (e) is 72 ± 12 events. 
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Figure 13: (from ref. 8) CDF measurement of the average B meson differential 
cross section based on the signals in the previous Figure. The branching ratio uncer- 
tainty, which contributes 11.9% to the cross section uncertainty, is shown separately. 
The dashed curves indicate the uncertainty in the theoretical prediction due to the 
uncertainties in m/,, momentum fragmentation and renormalization scale. The theo- 
retical prediction assumes f B + — f B o = .375 for the hadronization fractions of b (or 
b) quarks into the mesons indicated. The inset plot shows (data-QCD)/QCD on a 

linear scale. 
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Figure 14: (from ref. 9) The differential cross section for 8(f) from the CDF exper- 
iment. The angle 5<j> is measured in the transverse plane between the directions of 
the muon from a 6 quark and the jet from the b quark. The b jet must have a vertex 
tag, Et > 10 GeV and \r)\ < 1.5 . The muon must satisfy pr > 9 GeV and |?7| < 0.6. 
The data points have a common systematic error of ±9.5%. The uncertainty in 
the theoretical curves comes from the muonic branching fraction and fragmentation 

model. 
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with vertex detector hits. 
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Figure 15: (from ref. 9) Invariant mass for At — > J/tpA, J/tp — » /J, fi~ , A — > p7r~. 
The vertical scale is number of events per 10 MeV. The cuts are discussed in the 
text. The plot on the left is for muons with vertex detector hits; that on the right 
for muons without vertex hits. 

The cross section for Af, can be measured, under certain assumptions, by 
comparing the measured signal for Af, — > J/tpA with the signal for B d — > 
J/ipK®. Both signals have similar topologies in the detector. From the signal 
sizes and acceptances it is found that, 

a ■ Br(A b J/tpA) 



0.31 ±0.15 ±0.06. 



a ■ Br{B d -> J/«) 

a(A 6 )/a(B d ) = 0.1/0.375 



Assuming that, 
and using, 

Br(B d -» J/V^°) = 3.7 x 10~ 4 , 
which is consistent with the measured value, gives, 

Br(A b -> J/ipA) = 4.3 ± 2.1 ± 0.8 x 10~ 4 . 



(7) 

(8) 
(9) 

(10) 



To the extent that we consider the later Br to be reasonable, the assumed 
ratio of cross sections is also reasonable and hence provides a rough estimate 
of the Af, cross section. 
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5 Quarkonia 



5. 1 Theory 

The 'quarkonia' are bound states of a heavy quark and its anti- quark. Char- 
monium is a cc state and bottomonium is a bb state. The states are further 
characterized by their spin, orbital angular momentum and radial quantum 
numbers. The L — charmonium states are the J/tp and tp family; the L = 
bottomonium states are the T family. These L = states have quantum 
numbers J PC = 1 . 

The theory of the hadroprpduction of quarkonia is discussed in recent pa- 
pers by Braaten and Fleming Ed and Cho and LeibovichE3. The odd charge 
conjugation quantum number of the 1 states leads to a fundamental restric- 
tion on their production in gluon collisions. The simple process gg — > QQ(l ) 
is not allowed since gluons have C = — 1. However, a similar production pro- 
cess is allowed in which a gluon is radiated from a final state quark line, as 
shown in Figure |l6|(a) . This diagram is of order a| . If the intermediate state 
is a color singlet 1, this process is called direct color singlet production. But, 
the intermediate state need not be restricted to a color singlet. 

The intermediate state may be a either a color singlet or octet 8. Further- 
more, the intermediate state can be have the full range of quantum numbers 
2S+1 Lj. Thus, the general process, 

gg^QQ( 2S+1 L,i)^g + QQ{l—) (11) 

is possible for the production of the J/tp and T families. Examples of the 
intermediate state are 3 Si, 1 Sq, and 3 Pj. An example of direct color octet 
production is shown in the diagram of Figure |l6|(b). 

A third important process, also of order a|, involves the direct production 
of the x states which are characterized by L — 1 and the quantum numbers 
PC = ++. Since C = +, x can be produced directly from a two gluon collision, 
without the emission of a final state gluon, as shown in Figure |T~6|(c) . 

The xci {J — 1) an d Xci {J — 2) states have large branching fractions 
into J/tp + 7. Thus, these states can be studied experimentally by detecting 
the J/tp alone or both final state particles. However, the xco (J — 0) has a 
small branching fraction into J/ijj + 7 and hence cannot be easily studied. 

From inspection of the simple production diagrams considered above, it is 
evident that x production has the most straightforward theoretical treatment. 
There is no soft gluon emission required in the fundamental production process, 
and furthermore, there is no fragmentation process as for a state with open 
charm or beauty. Direct J/ip or T production is complicated by the emission 
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Figure 16: Feynman diagrams illustrating two gluon production of J = 1~~ 
quarkonia. (a) Direct production via a color singlet intermediate state with gluon 
emission, (b) Direct production via a color octet intermediate state with gluon 
emission, (c) Direct production of % and subsequent decay to 7 + 1~~ quarkonium. 
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of the soft gluon required by conservation of charge conjugation parity of the 
strong interaction. 



5.2 Charmonium 

The J I ip was co-discovered in a hadronic fixed target experiment and has since 
been studied in that way with high statistics. A number of results have been 
Reported recently, from Fermilab fixed target experiments using high energy tt~ 
c3 and proton I 14 " 15 ! beams. As an illustration, we consider the measurement 
of J '/ip and tp(2S) production from Fermilab experiment E789 which used the 
extracted 800 GeV proton beam. 



The invariant mass plot from E789 is shown in Figure 17. Of par- 

ticular note is the large number of signal events allowing for detailed studies 
of the Xf and Pt dependence of the production cross section. 

The differential cross section in Pt for J / ip is shown in Figure Also 
shown are predictions from a computer program of Mangano based on leading 
order QCD. Shown separately are the contributions to the cross section from 
direct production and from the decay of xc states. The sum of the contribu- 
tions seems to fit the data well. But actually the curves shown are 7 times the 
prediction! Thus, the theory predicts the shape reasonably well but severely 
underestimates the absolute scale of the cross section. This situation is even 
more dramatic for ij)(2S) production. 

The E789 measurement of da/dP T for tp(2S) is shown in Figure [n]. A 
multiplicative factor of 25 is required for the leading order prediction to match 
the scale of the data. Like the J/tp, the shape of the prediction for ip(2S) is in 
reasonable agreement with the data. 

As mentioned in the theory section, measurements of \c production are 
important because of the firmer theoretical understanding of x versus vector 
meson quarkonium production. Unfortunately, the measurement of the photon 
is typically inefficient in detecting the decay xc - > + 7- A review of fixed 
target data on xc production is given by McManus One simple prediction 
is that the ratio of cross sections for xci and XC2 is given by spin counting. 
Thus, one expects, 

a(xci) _ 2(1) + 1 



a( X C2) 2(2) + 1 



= 0.6 (12) 



Data for this ratio from various experiments are plotted in Figure g0[ Within 
the large statistical errors, the data are indeed consistent with the simple spin 
counting prediction. This result is consistent with the picture that leading 
order QCD processes are sufficient to explain x production. 
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Figure 17: (from ref. 15) Invariant mass of pairs from Fermilab E789. The 

dashed curves show the fits to the combinatoric background. The name tp' is in 

common usage for the ip(2S) 
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Figure 18: (from ref. 15) Measurement from Fermilab E789 of the differential cross 
section for J/tp production versus the transverse momentum of the J/tp. The 20% 
systematic normalization uncertainty is not shown. The dashed curves are 7 times the 
leading-order predictions for the inclusive Pt distributions of J/tp mesons originating 
from various quarkonium states. The solid curve is the sum of the quarkonium 
contributions. (The contribution from B — > J/tp + X decays is negligible.) 
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Figure 19: (from rof. 15) E789 measurement of da/dP? for pN -» V(2S) + X 
versus Pt of the ip(2S). Also shown is the measurement of J/ip for comparison. The 
systematic normalization uncertainty of 30% for the tp(2S) is not shown. The solid 
curve is 25 times the leading order prediction for the ip(2S). 
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There have also been recent measurements of charmonium production from 
the Tevatron Collider. The CDF collaboration has reported preliminary mea- 
surements E_ZI of the differential cross sections for J/tp and ip(2S) production. 
The data for the cross section times branching fraction for pp — » J/ip + X, 
J /ip — -» versus Pt of the J/ip is shown in Figure The contribution 

of J/ip from B decay is eliminated using vertex information. The contribu- 
tion of J j ip from xc decays is measured separately by explicit reconstruction 
of the decay \c ~^ J/ip + 7- The measured cross section for prompt J/ijj 
production (including the small contribution from xc decays) is given by the 
open triangles. The perturbative QCD prediction, which is dominated by the 
contribution from xc production, underestimates the data by a factor between 
5 and 10 over the Pt range measured. Thus, the situation at Ecm = 1-8 TeV 
is remarkably similar to that of the fixed target data at Ecm = 0.04 TeV. 
Note also that the CDF measurement of J/ip from xc decay is in reasonable 
agreement with the QCD prediction. 

The CDF measurements for ip(2S) production are shown in Figure p2[ 
The prompt ip{2S) data points are the solid squares. The QCD prediction 
for prompt production is the area between the dash-dotted lines. The data 
and theory disagree in both shape and normalization. At the lowest Pt the 
prediction is at least 20 times too low. As Pt increases, the discrepancy gets 
even larger. These results, along with those at fixed target energies, indicate 
that non-perturbative processes (such as the soft gluon emission shown in 
Figure |l6|) play an important role in prompt J/ip and ip(2S) production. 

Explicit—reconstruction of xc has also been accomplished with the DO 
experiment 113 at the Fermilab Tevatron. The signal can be seen in Figure 23 



which shows the distribution of the difference in invariant mass between the 
reconstructed J/ip — > and the candidate xc ~ * A* M~7 combinations. 

Note that the resolution is not sufficient to resolve the J = 0, 1, 2 xc states as 
is possible in fixed target experiments. It is found in DO that the fraction f x 
of J/ip's originating from xc is, 

f x = 0.32 ±0.07 ±0.07. (13) 

That this ratio is far from one indicates that the bulk of J/ip production is 
direct and not from xc decay. This reinforces the importance of the non- 
perturbative process of direct production. 

5.3 T Production 

The T was discovered using fixed target hadroproduction and has now been 
well measured at the Fermilab Tevatron collider. An invariant mass spectrum 
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Figure 21: (from ref. 17) CDF measurement of J/i/> production. Data is shown for 
the cross section times branching fraction for pp — ► J/t/j + X, J/?/) — ► versus 
Pt of the J/V>- The contribution of J /if) from B decay is eliminated using vertex 
information. The contribution of J/ip from \c decays is measured separately by 
explicit reconstruction of the decay \c — > J/tp + 7- The measured cross section for 
prompt J/tp production (including the small contribution from \c decays) is given 
by the open triangles. The perturbative QCD prediction, which is dominated by the 
contribution from \c production, underestimates the data by a factor between 5 and 
10 over the Pt range measured. 
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Figure 22: (from ref. 17) CDF measurement of tp(2S) production. The prompt tp(2S) 
data points are the solid squares. The QCD prediction for prompt production is the 
area between the dash-dotted lines. The data and theory disagree in both shape and 

normalization. 
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for T — » from CDfH is shown in Figure |24|. The experimental resolution 

is good enough to see separate peaks for the T IS*, 25 and 35 radial states. The 
differential cross sections in Pt for these particles are shown in Figure [2^. Also 
shown are the predictions of the color singlet model of Cho and Leibovitch. The 
prediction includes direct production and feeddown from xb- The predictions 
fall below the data for the three particles shown. For the T(15), the shape 
of the prediction is not in agreement with the data. Thus, it appears that 
the theoretical treatment of the production of the J = 1 bottomonium states 
suffers the same difficulty as for the J — 1 charmonium states: nonperturbative 
effects arc likely important. 

6 Concluding Remarks 

In this review, we've examined both open and closed charm and beauty states 
produced at CM energies spanning nearly two orders of magnitude (0.02 TeV 
to 1.8 TeV). Experimental accomplishments include measurements of fully re- 
constructed particles yielding normalized cross sections and their dependence 
on transverse and longitudinal momentum. For the open states, perturbative 
QCD calculations are successful in explaining the shapes of the differential 
cross sections but underestimate the scale of the cross sections by factors of 
2-3. Furthermore, the theoretical predictions are too sensitive to variations in 
the underlying parameters of the theory, such as the scales of renormalization 
and factorization. This leads to theoretical uncertainties in the cross section 
scale of order of the difference between the measurement and prediction. 

The picture emerging for quarkonia production is that for the J/ip and T 
families the production mechanism is fundamentally non-perturbative. This is 
evidenced by the large discrepancy between the measured cross sections and 
those predicted using QCD perturbation theory. However, measurements of \C 
production seem to be reasonably explained as due to perturbative processes. 

Future experiments promise progress in several areas. The Tevatron col- 
lider experiments will continue to study production of open beauty, and both 
charm and beauty quarkonia. LHC experiments will open a new energy regime. 
A new fixed target experiment, SELEX, is underway at Fermilab using X and 
proton beams. The COMPASS experiment is planned for the future CERN 
fixed target program using a hyperon beam. The HERA-B experiment is under 
construction at DESY and will study fixed target pN production. 

The challenge to theorists is to tame the theoretical errors in perturbative 
QCD calculations with new approaches to factorization and renormalization. 
Models that can treat the fundamentally non-perturbative processes in quarko- 
nia production and in fragmentation of heavy quarks would be an important 
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Figure 24: (from ref. 19) CDF measurement of the invariant mass distribution of 
muon pairs in the T mass region for \y\ < 0.4. There are also cuts on the Pt of the 
muons. The histogram corresponds to the data and the solid curve is a fit to three 
Gaussians plus a quadratic background. 
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Figure 25: (from ref. 19) CDF data for T production compared to the predictions 
of the QCD color singlet model of Cho and Leibovitch. Plotted is the product of the 
cross section times branching fraction (d 2 a/dPTdy) y= o x Br versus Pt for T. The 
predictions (solid lines) include direct production and feeddown from \b- In c), the 
dashed line corresponds to direct production only, whereas the solid line also includes 
the contribution from the decay of the unobserved Xb(3P) state. 
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tool to assess whether "new" physics is required. 

The hadroproduction of charm and beauty particles serves as one of the 
testing grounds for QCD. Hadronic collisions also provide a copious source of 
heavy quarks, allowing us to study their spectroscopy and decay with great 
precision. However, we must understand the production characteristics well 
to achieve precision in the decay measurements. These incentives spark our 
interest in future work in this field. 
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